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摘 要 : 基于 MODIS-EVI 数据 ,采用 一 元 线性 回归 经 验 正 交 函数 (EOF) 、 变 异 系数 ,从 时 间 和 空间 2 个 维度 分 析 
古 尔 班 通 古 特 沙漠 增强 型 植被 指数 (EVI) 的 时 空 变化 特征 。 结 果 表 明 :2000 一 2018 年 古 尔 班 通 古 特 沙漠 整体 EVI 


的 年 际 变化 呈 显 著 增加 趋势 ,增长 速率 为 0.016 0. (10a) ^ (P <0.01) ,其 中 国定 沙丘 、 半 固定 沙丘 和 流动 沙丘 中 
的 EVI 年 际 变化 也 呈 显 著 增加 趋势 ;生长 季 内 沙漠 整体 及 不 同类 型 EVI 的 变化 趋势 大 致 相同 ,从 3 


omn ,并 在 6 一 7 


和 开始 ,EVI 3 


月 达到 一 年 中 的 最 高 值 ; 古 尔 班 通 古 特 沙漠 大 部 分 区 域 的 EVI 呈 上 升 趋势 ,上 升 趋势 明显 的 区 


域 主要 集中 在 沙漠 南 缘 以 及 沙漠 西 缘 的 开荒 区 域 , 沙漠 腹地 EVI 上 升 趋势 较 小 ,EVI 降低 的 区 域 主要 分 布 在 沙漠 
北部 ; 古 尔 班 通 古 特 沙 漠 EVI 波动 较 大 的 区 域 在 沙漠 西 缘 与 南 缘 , 其 中 波动 最 大 的 区 域 在 沙漠 西 缘 的 开荒 区 域 , 沙 


漠 腹 地 EVI 波动 较 小 。EVI 能 够 较 好 地 反映 以 固定 和 半 固 


域 沙 漠 地 貌 类 型 的 空间 划分 具有 重要 的 参考 意义 。 


定 为 主 的 沙漠 区 植被 覆盖 变化 ,其 反映 的 植被 状况 对 区 


关键 词 : 增强 型 植被 指数 (EVI1) ; 时 空 变化 ; 经 验 正 交 函数 (EOF) ; 古 尔 班 通 古 特 沙漠 


我 国 是 世界 上 干旱 地 区 和 沙漠 分 布 最 多 的 国家 
之 一 。 沙 漠 地 区 的 开发 利用 、 建 设 与 保护 对 区 域 
环境 与 经 济 的 发 展 具 有 重要 意义 。 

古 尔 班 通 古 特 沙漠 是 我 国 最 大 的 固定 半 固 定 沙 
汉中 ,具有 独特 的 沙漠 地 貌 和 植被 ,前 人 研究 主要 
集中 在 植物 群落 、 区 系 及 其 分 布 等 方面 *“ , 且 大 部 
分 是 采用 传统 的 野外 调查 方法 ,范围 有 限 、 时 间 序 列 
短 。 随 着 各 感 技术 的 发 展 ,采用 遥感 手段 监测 植被 
生长 状况 成 为 目前 研究 的 主流 方式 ,遥感 技术 有 助 
于 监测 生境 季节 变化 和 年 际 变化 ,可 以 进行 时 间 序 
列 分 析 。 一 些 学 者 利用 长 时 间 序 列 遥 感 数据 ,对 古 
尔 班 通 古 特 沙漠 植被 的 时 空 变 化 进行 了 研究 。 李 杨 
等 “分 析 了 1981—2003 年 古 尔 班 通 古 特 沙漠 不 同 
分 区 ,基于 归 一 化 植被 指数 (normalized difference 
vegetation index, NDVI) 时空 变化 特征 ; 蒋 超 亮 等 
根据 2000 一 2016 年 NDVI 数据 提取 了 该 沙漠 早生 
植物 的 时 空格 局 ; 段 呈 等 ”对 1985 一 2014 年 的 短命 
植物 NDVI 的 时 空 变化 进行 了 研究 。 以 往 的 研究 大 
都 基于 NDVI 进行 分 析 , 考 虑 到 使 用 传统 的 归 一 化 
植被 指数 存在 许多 不 足 和 缺陷 ,而 增强 型 植被 指数 
(enhanced vegetation index, EVI) 能 够 减少 来 自 大 气 
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和 土壤 噪音 的 影响 17 ,并 且 EVI 时 间 序 列 相 较 
于 NDVI 时 间 序 列 季节 性 变化 更 明显 ,时 间 序 列 曲 
Ag, 

沙漠 地 貌 复 杂 而 多 样 , 存 在 各 种 基本 地 貌 类 型 、 
成 因 类 型 ,形态 类 型 ,其 表层 覆盖 有 不 同类 型 和 盖 度 
的 植被 ,并 随 着 自然 及 人 为 的 影响 发 生 着 变化 ,通过 
对 沙漠 植被 的 监测 有 益 于 多 角度 地 认 知 沙漠 。 本 研 
究 选用 增强 型 植被 指数 (EVI) ,采用 一 元 线性 回归 
分 析 、 经 验 正 交 函 数 (empirical orthogonal function, 
EOF) 变异 系数 ,对 2000—2018 年 古 尔 班 通 古 特 沙 
漠 EVI 的 时 空 变化 特征 进行 研究 ,这 对 区 域 沙 漠 地 
貌 类 型 的 认识 具有 重要 的 参考 意义 ,可 为 沙漠 分 类 、 
制图 和 数据 库 的 构建 提供 基础 数据 ,并 且 对 于 沙漠 
的 动态 监测 .保护 及 利用 也 具有 十 分 重要 的 意义 。 


1 研究 区 与 数据 来 源 
1.1 研究 区 概况 

古 尔 班 通 古 特 沙漠 位 于 新 疆 北部 准噶尔 盆地 腹 
地 ,总 面积 约 为 4.73 x 10* km ,是 我 国 第 二 大 沙漠 
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(图 1) 。 沙 漠 南 缘 地 接 天 山北 坡 冲 洪 积 平原 区 的 绿 
洲 荡 漠 带 ,多 靠近 天 山北 芳 诸 多 河流 的 下 游 尾 闻 ,水 
土 条 件 较 好 ,是 新 疆土 地 开发 程度 较 高 的 区 域 ; 沙 漠 
西 缘 是 玛 纳 斯 河 下 游 、. 玛 纳 斯 湖 盆 区 , 玛 纳 斯 河 、 塔 
西河 - 呼 图 壁 河 老 河道 穿越 沙漠 西南 部 ,也 是 农垦 
快速 发 展 地 区 ;沙漠 北部 是 乌 伦 古 河南 部 荒地 ,主要 
为 干燥 剥蚀 基底 平原 ,无 天 然 河流 ,气候 环境 较 差 ; 
沙漠 东部 是 准噶尔 东部 的 低 山 .丘陵 及 平原 区 ,干燥 
剥蚀 和 风蚀 强烈 。 

古 尔 班 通 古 特 沙漠 区 年 降水 量 可 达 70 ~ 150 
mm ,7 一 8 月 降水 较 多 ,冬季 有 稳定 积 雪 , 积 雪 的 消 
融 下 渗 为 春季 短命 植物 的 生长 提供 了 有 利 的 条 
件 小 。 固 定 、 半 固定 沙丘 是 古 尔 班 通 古 特 沙漠 的 主 
要 地 貌 类 型 ,固定 沙丘 主要 分 布 在 沙漠 西部 .南部 边 
缘 以 及 生态 恢复 的 区 域 , 植 被 覆盖 度 大 于 3096 ; 半 
国定 沙丘 植被 覆盖 度 在 10% ~ 30% ,主要 分 布 于 沙 
漠 腹 地 ;流动 沙丘 植被 覆盖 度 小 于 10% , 主要 分 布 
于 沙漠 的 东部 边缘 ,广泛 发 育 在 沙丘 顶部 ,多 为 流动 
EPES, 

古 尔 班 通 古 特 沙漠 内 部 植物 长 势 较 好 ,植被 组 
成 相对 丰富 ,其 植物 区 系 和 植被 性 质 表 现 出 蒙古 蕊 
壁 荡 漠 与 中 亚 西 部 荒漠 之 间 的 过 渡 特 征 , 具 有 世界 
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温带 荒漠 中 最 为 丰富 的 植物 物种 与 基因 资源 。 
在 春季 ,利用 降水 和 积 雪 融 水 发 育 的 短命 植物 .类 短 
命 植物 是 古 尔 班 通 古 特 沙 漠 植 物 区 系 特有 的 现象 。 

古 尔 班 通 十 特 沙漠 的 半 固 定 沙 丘 主 要 分 布 着 白 
RIR . 贰 属 等 荒漠 植被 ;固定 沙 芍 沙 丘 上 主要 分 布 着 
HIRI ER ,麻黄 ,短命 植物 ;沙漠 东 缘 的 半 流 动 沙 
丘 上 主要 为 沙 抛 枣 从 和 巨 穗 滨 麦 从 ; 玛 纳 斯 河 旁 则 
为 胡杨 林 、 盐 化 草 和 甸 、 河 漫 滩 草 甸 "” 。 近 年 来 河流 
下 游 及 沙漠 边缘 开展 严重 ,农作物 及 人 工 植被 扩展 
迅速 。 
1.2 数据 来 源 

本 研究 采用 的 MODIS 产品 来 源 于 美国 国家 航 
空 航天 局 (http ://ladsweb. nascom. nasa. gov ) 提供 的 
合成 植被 指数 产品 MOD13Q1, 空间 分 辨 率 为 250 
m x250 m, 时 间 分 辨 率 为 16 d, 时 空 跨度 为 2000 年 
3 月 至 2018 年 11 月 。 古 尔 班 通 古 特 沙漠 研究 区 共 
2 景 影像 ,其 过 境 轨道 号 分 别 为 h23v04 和 h24v04。 
MODIS 琐 感 数据 预 处 理 使 用 美国 国家 航空 航天 局 
(National Aeronautics and Space Administration, 
NASA ) 免 费 提 供 的 MRT(MODIS reprojection tool) 软 
件 , 预 处 理 包括 影像 拼接 ,波段 提取 、 重 投影 .格式 转 
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注 : 右 上 角 图 基于 新 疆 维吾尔 自治 区 测绘 地 理 信息 局 标准 地 图 服务 网 站 下 载 的 审 图 号 为 新 S(2016)148 号 的 标准 地 图 制作 , 底 图 无 修改 。 
图 1 上 古 尔 班 通 古 特 沙漠 类 型 及 其 在 新 疆 的 位 置 


Fig.1 The types and location of the Gurbantunggut Desert in Xinjiang 
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研究 区 沙漠 类 型 数据 来 源 于 中 国 科学 院 资源 环 
境 与 地 理 研 究 所 .中国 科 学 院 新 疆 生 态 与 地 理 研究 
所 等 单位 编制 的 《中 华人 民 共 和 国 地 貌 图 集 (1: 100 
万 )》, 从 中 裁 切 出 古 尔 班 通 古 特 沙漠 区 域 的 空间 数 
据 " (图 1)。 


2 研究 方法 


2.1 Savitzky-Golay 滤波 

EVI 时 序数 据 集 须 在 应 用 前 正确 有 效 地 去 除 噪 
声 , 重 建 EVI 时 序数 据 集 "1。 采 用 Savitzky-Golay 
滤波 进行 EVI 时 间 序 列 重 构 ,能够 很 好 地 对 EVI 时 
序曲 线 进行 平滑 ,反映 植物 生长 变化 趋势 ,并 且 具 有 
实现 简单 . 先 验 知识 少 等 优点 7。 

S-G 滤波 是 Savitzky 和 Golay 提出 的 基于 多 项 
式 最 小 二 乘法 拟 合 的 一 种 数据 平滑 方法 ,可 以 将 
其 理解 为 一 种 加 权 滑 动 平均 滤波 ,其 权重 取决 于 一 
个 滤波 窗口 范围 内 做 多 项 式 最 小 二 乘法 拟 合 的 多 项 
式 次 数 。S-G 滤波 过 程 可 以 用 公式 (1 HRI 


M 
Ee (1) 

式 中 :六 为 合成 序列 数据 即 EVI 拟 合 值 ;m 为 窗口 
宽度 ,定义 为 平均 窗口 大 小 的 1/2;C; 为 多 项 式 拟 合 
的 系数 ,表示 从 滤波 器 首部 开始 第 i 个 EVI 值 的 权 
重 ;7,,; 代 表 原 始 序 列 数 据 即 EVI 原始 值 ;N 为 滑动 
窗口 所 包括 的 数据 点 ( 即 滤波 器 的 长 度 ) ,等 于 滑动 
数组 的 宽度 (2m +1) 。 在 本 研究 中 ,将 S-G 滤波 的 
和 窗口 设置 为 4, 过 代 次 数 设 置 为 3, 对 2000 一 2018 年 
古 尔 班 通 古 特 沙漠 EVI 时 间 序 列 数据 进行 去 噪声 
拟 合 ,以 重 构 EVI 时 间 序 列 数据 集 。 
2.2 线性 回归 分 析 

采用 一 元 线性 回归 方法 ,分 析 古 尔 班 通 古 特 沙 
漠 各 年 份 年 最 大 化 EVI 的 变化 趋势 ,以 单个 像 元 时 
间 特 征 反应 整个 区 域 趋势 变化 。! 是 利用 最 小 二 乘 
法 公式 [公式 (2) ] ,基于 每 个 像 元 计算 的 年 最 大 化 
EVI 变化 趋势 。 


žen- a e) 
»u i ( ZEVI, ) [N 
AP N 为 研究 范围 的 总 年 数 ,本 文中 W 取 19;i 为 


人 研究 期 中 的 第 i 年 ;EVI; 为 第 i 年 的 年 最 大 合成 
EVI。 某 像 元 点 的 趋势 线 是 这 个 像 点 19 a 的 年 最 大 


(2) 


化 EVI 值 用 一 元 线性 回归 模拟 出 来 的 一 个 总 的 变 
化 趋势 ,5 为 这 条 趋势 线 的 斜率 ,斜率 为 负 表 示 研 究 
区 内 EVI 变化 趋势 为 下 降 , 反 之 , 则 表示 研究 区 内 
EVI 变化 趋势 为 上 升 。! > 0 时 ,数值 越 大 ,表示 EVI 
增加 趋势 越 明 显 ; <0 时 ,数值 越 小 ,表明 EVI 降低 
趋势 越 明 显 。 变 化 趋势 的 显著 性 采用 /上 检验 。 
2.3 经 验 正 交 函数 分 解 

经 验 正 交 函 数 (EOF) 分 解 是 一 种 在 地 球 科学 领 
域 得 到 广泛 应 用 的 场 分 析 方 法 ,其 原理 是 将 时 空 要 
素 场 分 解 为 若干 空间 的 基本 模 态 和 时 间 系 数 序 列 的 
线性 组 合 , 进 而 客观 定量 分 析 要 素 场 的 空间 结构 和 
时 间 变 化 宇 -六 。 在 本 研究 中 用 以 分 析 古 尔 班 通 十 
特 沙漠 EVI 的 时 空 变化 特征 。 设 研究 区 共有 m 个 
观测 点 ,每 个 点 有 nn 次 观测 值 , 则 观测 数据 可 写成 公 
式 (3) 的 矩阵 形式 。 通 过 EOF 展开 ,将 公式 (3) 分 
解 成 正 交 的 空间 和 抢 阵 了 与 正 交 的 时 间 和 抑 阵 了 乘积 
之 和 [公式 (4) ] 。 


X= (X,,) = (XiX; Xn) (3) 
式 中 XX, X, jS 为 维 数 是 m 的 列 向 量 。 
X-VT (4) 


空间 矩阵 了 可 由 XX” 的 特征 向 量 得 出 [公式 
(5) ] ,时 间 和 矩阵 由 公式 (6) 可 计算 。 
C - XX! = VIT V = VAV" (5) 
T-V'X (6) 
式 中 :上 和 角 标 7” 表示 和 矩阵 转 置 ;4 为 XX" 矩阵 的 特 
征 值 构成 的 对 角 阵 。 
空间 模 态 的 贡献 率 R, 可 由 公式 (7) 计算 得 出 。 


2 (7) 
à 
式 中 :RR, 为 第 个 模 态 的 方差 贡献 率 ;A A XX" 的 
特征 值 ;A; 为 第 i 个 模 态 的 方差 贡献 ;m 为 列 向 量 的 
维 数 。 

在 进行 结果 分 析 前 还 需要 检验 分 解 出 来 的 空间 
模 态 是 否 有 物理 意义 ,本 文采 用 特征 值 误 差 范 围 的 
方法 来 进行 检验 。 空 间 模 态 的 贡献 率 越 大 , 越 能 
反映 变量 场 的 空间 分 布 特征 ,空间 模 态 数值 的 绝对 
值 越 大 ,说 明 其 随 着 时 间 的 变化 幅度 越 大 。 时 间 系 
数 是 空间 模 态 的 权重 ,描述 空间 模 态 的 时 间 变 化 信 
息 ,数值 越 大 ,表明 对 应 年 份 的 空间 分 布 特征 与 该 空 
间 模 态 越 相似 。 
2.4 波动 特征 分 析 

变异 系数 (coefficient of variance ,下 文 以 Cv 表 


,bz1,2, 5, p(p«m) 
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示 ) 是 衡量 资料 中 各 观测 值 变 异 程度 的 统计 量 ,是 
标准 差 与 均值 的 比率 ,可 消除 单位 或 平均 数 对 不 同 
2 个 或 多 个 资料 变异 程度 比较 的 影响 ,以 反映 单位 
均值 的 离散 程度 人 “”?。 本 研究 通过 计算 变异 系数 
[公式 (8)] 来 分 析 2000—2018 年 古 尔 班 通 古 特 沙 
iE EVI 的 波动 特征 。 

Cv = 二 (8) 

H 

xp io 为 标准 差 ;y 为 平均 值 。 采 用 变异 系数 作为 
波动 性 评价 指标 ,分 析 古 尔 班 通 古 特 沙漠 EVI 的 波 
动 规律 ,Cv 值 越 大 ,表示 EVI 值 的 分 布 越 离散 ,波动 
较 大 ;反之 ,表示 EVI 值 分 布 越 集中 ,波动 较 小 。 


3 结果 与 分 析 


3.1 EVI 的 时 空 分 布 

基于 2000 一 2018 年 EVI 值 的 年 最 大 值 合成 法 ， 
计算 得 到 其 19 a 平均 值 的 空间 分 布 图 ,由 此 分 析 得 
到 : 古 尔 班 通 古 特 沙 漠 19 a 来 的 EVI 年 最 大 均值 在 
0.04 ~0.87 之 间 ,EVI 值 较 高 的 地 区 主要 分 布 在 沙 
漠 的 西南 部 和 南 缘 的 沙漠 - 绿洲 过 渡 带 ,沙漠 腹地 
的 EVI 值 较 沙 漠 的 西 东 两 侧 EVI 值 高 (图 2)。 


EVI 
Z 0.87 
Ü - 
0.04 
Z 
g- 
< 
z 0 30 60 120 km 
E 
+ 


85°E 86°E 87°E 88°E 89*E 90E — 91^E 

到 2 2000—2018 年 古 尔 班 通 古 特 沙漠 EVI 空间 分 布 

Fig.2 The spatial distribution of EVI in the Gurbantunggut 
Desert from 2000 to 2018 


3.2. EVI 的 年 际 变化 

基于 年 最 大 化 合成 EVI 数据 ,采用 一 元 线性 回 
归 分 析 方法 对 2000—2018 年 十 尔 班 通 十 特 沙漠 及 
其 不 同 沙漠 类 型 上 EVI 的 年 际 变化 特征 进行 分 析 。 
结果 表明 (图 3) : 古 尔 班 通 古 特 沙漠 不 同 沙漠 类 型 


一 古 尔 班 通 古 特 沙漠 y-0.0016x-3.097 1, P « 0.01 


-- 半 固 定 沙丘 y=0.001 6x—-3.100 9, P< 0.01 

oio- = 固定 沙丘 y-0.001 6x—3.009 7, P < 0.01 

; 一 流动 沙丘 y=0.000 8x—-1.415 5, P< 0.01 
0.145 + 
= 0.120 上 

0.095 上 “ 

[5] © © Q =] ml Dm — = 一 

[-] [=] =] =] © > © © 

N N N N N N N N N N 


3 2000—2018 年 古 尔 班 通 古 特 沙漠 EVI 的 年 际 变化 
Fig.3 Annual change of EVI in the Curbantunggut Desert 
from 2000 to 2018 


0.105 0 出 现在 2001 年 ,最 高 值 0. 163 3 出 现在 
2016 年 ,上 升幅 度 为 55. 5296 ,表明 近 19 a 来 , 古 尔 
班 通 古 特 沙漠 EVI 呈 增 长 态势 ,增长 速率 较 快 。 半 
固定 沙丘 和 固定 沙丘 EVI 的 年 际 变化 趋势 与 整个 
古 尔 班 通 古 特 沙漠 EVI 的 变化 趋势 相似 ,其 中 半 固 
定 沙丘 EVI 以 0.016 1 . (10a) ”的 变化 速率 呈 显 
著 增 长 (P <0.01) ,EVI 变化 在 0.104 8 ~0. 163 5 
之 间 ,2016 年 的 EVI 值 0.163 5 为 整个 研究 时 间 段 
的 最 大 值 ; 古 尔 班 通 古 特 沙漠 固定 沙丘 EVI 的 变化 
速率 为 0.015 6 - (10a) -!(P <0.01) , 略 小 于 半 轩 
定 沙丘 EVI 的 变化 速率 , EVI 变化 在 0. 099 1 ~ 
0.155 0 之 间 , EVI 最 高 值 0. 155 0 出 现在 2016 年 。 
流动 沙丘 的 EVI 变化 速率 最 低 , 为 0. 007 6 - 
(10a) "(P «0.01) ,EVI 变化 在 0.093 1 ~0.1194 
之 间 ,EVI 最 大 值 0.119 4 出 现在 2016 年 ,EVI 最 小 
值 0.093 1 出 现在 2001 年 。 
3.3 ”EVI 的 生长 季 变 化 

基于 整个 古 尔 班 通 古 特 沙 漠 EVI 月 均值 ,分 析 
其 植物 生长 季 (3 一 10 月 ) EVI 的 年 内 变化 (图 A), 
结果 表明 :生长 季 内 十 尔 班 通 古 特 沙 漠 不 同 沙漠 类 
型 EVI 值 的 变化 趋势 大 致 相同 ,从 3 月 开始 ,EVI 值 
逐渐 增加 ,并 在 6 一 7 月 达到 一 年 中 的 最 大 值 。 各 类 
型 中 半 固 定 沙丘 EVI 值 最 大 ,固定 沙丘 EVI 值 次 


EVI 均 呈 现 显 著 增加 的 趋势 (P<0.01), 且 古 尔 班 
通 古 特 沙漠 固定 沙丘 EVI 值 > 半 固 定 沙 丘 EVI 
值 > 流动 沙丘 EVI 值 。 从 整个 沙漠 来 看 , 古 尔 班 通 
古 特 沙漠 EVI 的 增长 速率 为 0.016 0 + (102) ^, 
EVI 变 化 在 0.1050 ~0.1633 之 间 ,EVI 最 低 值 


之 ,流动 沙丘 FEVI 值 最 小 。 古 尔 班 通 上 古 特 沙漠 的 
EVI 值 从 3 月 的 0.097 8 快速 增加 ,7 月 达到 峰值 ， 
峰值 处 EVI 值 为 0. 110 4,8 月 开始 EVI 有 所 下 降 ， 
10 月 降 至 最 小 , EVI 最 小 值 为 0.051 8。 半 固定 类 
型 沙丘 EVI 的 年 内 变化 特征 与 整个 古 尔 班 通 古 特 
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Fig.4 The change of EVI in different desert types of the 


Gurbantunggut Desert in growth season from 2000 to 2018 


沙漠 的 年 内 变化 基本 一 致 ,EVI 最 大 值 0. 110 4 出 现 
在 7 月 ,EVI 最 小 值 0.052 0 出 现在 10 月 。 古 尔 班 通 
古 特 沙 漠 中 国定 沙丘 各 月 份 的 EVI 值 变化 在 
0.042 9 ~0.106 6 之 间 ,EVI 值 由 3 月 的 0.092 2 逐渐 
增 大 ,在 6 月 达到 峰值 ,峰值 处 EVI 值 为 0.106 6 , 随 
后 EVI 值 逐 渐 下 降 至 最 小 。 流 动 沙 丘 EVI 值 年 内 变 
化 较 小 ,变化 范围 在 0.055 1 ~0. 096 6 之 间 , 各 月 份 
的 EVI 值 均 小 于 固定 沙 乒 和 半 固 定 沙丘 的 EVI 值 。 
3.4 EVI 的 线性 趋势 分 析 与 EOF 分 析 
古 尔 班 通 古 特 沙漠 EVI 数据 的 EOF 分 析 , 第 一 
空间 模 态 的 方差 贡献 率 为 96. 86% ( 表 1) , 且 通 过 
寺 征 根 误差 范围 显著 性 的 检验 ,能 够 较 好 的 表示 
2000—2018 年 古 尔 班 通 古 特 沙 漠 EVI 变量 场 的 主 
要 空间 分 布 特征 。 在 第 一 空间 模 态 中 (图 5) , 古 尔 
班 通 古 特 沙 漠 的 数值 均 为 正 值 ,在 空间 上 具有 很 好 
的 一 致 性 ,并且 数 值 的 高 值 区 主要 分 布 在 古 尔 班 通 
古 特 沙漠 的 南 缘 与 中 部 区 域 , 结 合 其 时 间 系 数 可 以 
发 现 , 在 2000—2018 年 , 古 尔 班 通 古 特 沙漠 的 EVI 
整体 呈现 出 上 升 的 趋势 ,其 中 古 尔 班 通 古 特 沙 漠 的 


X1 基于 古 尔 班 通 古 特 沙漠 EVI 数据 的 经 验 正 交 分 解 前 
2 个 模 态 的 方差 贡献 率 
Tab.1 The variance contribution rates of the first two 
modes of EVI in the Gurbantunggut Desert from 


empirical orthogonal function 


模 态 方差 贡献 /% 累计 方差 贡献 /% 
1 96. 86 96. 86 
2 2.30 99.16 
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图 5 2000—2018 年 古 尔 班 通 古 特 沙 漠 EVI 数据 EOF 
分 解 的 第 一 空间 模 态 及 其 对 应 的 时 间 系 数 


Fig.5 The first spatial mode and its corresponding time 


coefficients obtained from empirical orthogonal function 


of EVI in the Gurbantunggut Desert from 2000 to 2018 


南 缘 与 中 部 区 域 EVI 的 上 升 趋势 最 为 明显 。 基 于 
像 元 的 线性 趋势 分 析 结 果 与 EOF 分 析 第 一 空间 模 
态 的 结果 基本 一 致 (图 6, 表 2) , 古 尔 班 通 古 特 沙 漠 
大 部 分 区 域 (96.46% ) 的 EVI RE EIR, A 
3.54% 区 域 的 EVI 呈 降低 趋势 ,有 43% 的 区 域 显 著 
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图 6 2000—2018 年 古 尔 班 通 古 特 沙漠 EVI 时空 变 化 


Fig.6 Spatial change of EVI in the Curbantunggut Desert from 2000 to 2018 
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表 2 上古 尔 班 通 古 特 沙漠 EVI 年 际 变化 趋势 统计 
Tab.2 Statistic results of EVI trend in the 


Gurbantunggut Desert / 96 
ET 1 所 占 面积 比例 ”所 占 面积 比例 
变化 趋势 ”所 占 面积 比例 (P «0.05) (P «0.1) 
EFt 96.46 28.91 43.00 
下 降 3.54 0.05 0. 06 


mm 0.55-1.56 
=E 0.33-0.55 
== 0.20-0.33 
= 0.13-0.20 
= 0.03-0.13 
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图 7 2000—2018 年 古 尔 班 通 古 特 沙漠 EVI 的 
变异 系数 分 布 
Fig.7 The distribution of variable coefficient of EVI in the 
Gurbantunggut Desert from 2000 to 2018 


上 升 (P<0.1) ,0.06% 的 区 域 显著 降低 (P <0.1)， 
其 中 EVI 上 升 趋势 明显 的 区 域 主要 集中 在 古 尔 班 
通 古 特 沙漠 的 南 缘 , 上 升 趋势 最 为 明显 的 区 域 在 沙 
漠 西 部 边缘 ,沙漠 腹地 EVI 上 升 趋势 不 明显 ,EVI 降 
低 的 区 域 主要 分 布 在 沙漠 北部 。 
3.5 EVI 的 波动 特征 分 析 

基于 2000—2018 年 古 尔 班 通 古 特 沙漠 年 最 大 
化 合成 EVI 数据 ,采用 变异 系数 ( Cv) 分 析 古 尔 班 通 
古 特 沙 漠 植 物 生长 状况 的 波动 程度 (图 7)。 结 果 表 
8j :2000—2018 年 古 尔 班 通 古 特 沙漠 EVI 的 变异 系 
数 主要 集中 在 0.03 ~1.56 之 间 ,根据 自然 断 点 法 将 
其 分 为 5 级 ,0.03 < Co <0. 13 的 区 域 占 古 尔 班 通 古 
特 沙漠 的 48. 8996 ;0. 13 < Cv «0. 20 的 区 域 占 古 尔 
班 通 古 特 沙漠 的 30. 34% ;0. 20 < Cv «0. 33 的 区 域 
占 古 尔 班 通 古 特 沙 漠 的 17. 97% ;0. 33 < Cv <0. 55 
的 区 域 占 古 尔 班 通 古 特 沙漠 的 1.65% ;0. 55 < Cv < 
1.56 的 区 域 占 古 尔 班 通 古 特 沙 漠 的 1.15% 。 古 尔 
班 通 古 特 沙 漠 EVI 波动 程度 较 大 的 区 域 集中 在 沙 
漠 南 缘 与 西 缘 ,其 中 波动 最 大 的 区 域 在 沙漠 西 缘 的 
开荒 区 域 , 沙 漠 腹 地 及 北部 EVI 波动 较 小 。 


4 讨论 


4.1 EVI 指数 的 区 域 应 用 
EVI 在 一 定 的 环境 背景 影响 下 ,对 植被 变化 的 


敏感 性 要 强 于 NDVI,EVI 更 能 反映 区 域内 植被 的 空 
间 差 异 。 从 本 研究 可 以 看 出 ,EVI 可 以 应 用 在 沙漠 
区 域 , 且 能 够 较 好 表征 研究 区 植被 覆盖 状况 的 变化 
趋势 。 从 2000—2018 年 古 尔 班 通 古 特 沙 漠 EVI 的 
年 际 变化 趋势 可 以 看 出 , 近 19 a 来 古 尔 班 通 古 特 沙 
漠 及 其 不 同 沙漠 类 型 植被 覆盖 状况 均 向 好 的 趋势 发 
E RSF H WEA BE, TRET 
NDVI\RSEI( 遥感 生态 指数 ) 等 数据 得 出 的 研究 结 
果 基 本 一 致 。 

研究 小 组 于 2018 年 5 一 9 月 对 古 尔 班 通 十 特 沙 
漠 区 进行 了 实地 调查 ,包括 定点 样 方 调查 .无 人 机 空 
中 俯 辐 估算 植被 盖 度 .线路 沿途 定点 植被 盖 度 记录 
等 ,调查 发 现 :其 各 地 的 植被 盖 度 状况 与 基于 EVI 
得 出 的 植被 覆盖 状况 的 空间 格局 相 一 致 。 
4.2 EVI 指数 与 沙漠 地 貌 类 型 

目前 ,沙漠 地 貌 类 型 的 分 类 多 基于 地 貌 成 因 类 
型 与 地 貌 形态 类 型 。 本 研究 采用 的 沙漠 类 型 空间 分 
布 数据 由 于 制图 尺度 .数据 裁 切 ,沙漠 地 貌 分 类 方法 
等 的 差异 ,没有 将 该 沙漠 南 缘 的 许多 固定 沙丘 区 域 
纳入 空间 分 析 ,致使 固定 沙丘 整体 的 EVI 值 偏 小 。 
通过 野外 实地 调查 可 以 发 现 :2000 一 2018 年 古 尔 班 
通 古 特 沙漠 EVI 空间 分 布 图 可 以 较 好 地 反映 沙漠 
不 同 地 理 区 位 的 植被 状况 。 因 此 ,在 以 后 的 沙漠 分 
类 及 制图 中 ,可 以 将 基于 EVI 等 植被 指数 和 地 面 植 
被 调查 数据 得 出 的 植被 覆盖 状况 作为 划分 沙漠 地 貌 
类 型 的 重要 参数 。 
4.3 植被 覆盖 变化 分 析 

2000—2018 年 古 尔 班 通 古 特 沙漠 大 部 分 区 域 
的 EVI 都 旦 上 升 的 趋势 , 极 少 部 分 区 域 的 EVI 呈 降 
低 趋势 ;各 沙漠 类 型 区 的 EVI 均 呈现 显著 增加 趋势 
(P<0.01)。 气 温和 降水 是 影响 地 表 植 被 覆盖 状况 
的 2 个 关键 因素 ,水 分 条 件 尤为 重要 。 沙 漠 南 缘 和 
沙漠 东部 靠近 准噶尔 盆地 东部 低 山 的 区 域 降水 量 较 
多 ,沙漠 最 北部 风 多 温度 低 。 基 于 EVI 与 降水 的 偏 
相关 分 析 , 古 尔 班 通 古 特 沙漠 大 部 分 区 域 的 植被 覆 
盖 度 与 降水 呈正 相关 ,通过 对 准噶尔 区 域 气象 资料 
的 分 析 以 及 有 关 该 区 域 降水 气温 .干旱 等 研究 文献 
表明 :20 世纪 80 年 代 中 期 以 来 降水 量 普遍 增多 ,年 
降水 量 及 年 降水 日 数 均 呈 现 增加 趋势 , 区域 气 候 有 
向 暖 湿 发 展 的 趋势 ,这 为 古 尔 班 通 古 特 沙 漠 内 植物 
的 生长 提供 了 有 利 条 件 ” T, EF EVI 与 温度 
的 偏 相关 分 析 也 表明 , 沙漠 植被 覆盖 度 变化 与 气温 
存在 不 同 程度 的 正 ` 负 相关 关系 , 既 有 增长 也 有 抑制 


作用 。 沙 漠 区 域 植被 覆盖 的 变化 与 降水 -温度 协同 
的 量化 关系 有 待 进一步 研究 。 

古 尔 班 通 古 特 沙 漠 西 南部 与 南 缘 位 于 沙漠 — 绿 
洲 过 渡 带 ,地 下 水 位 较 高 ,水 热 条 件 相 对 较 好 ,植被 
覆盖 状况 变 好 的 趋势 最 明显 525) 。 研 究 区 植被 覆盖 
状况 变 好 的 原因 主要 有 3 个 方面 :其 一 , 古 尔 班 通 古 
特 沙漠 南 缘 及 西南 部 地 区 靠近 天 山北 坡 经 济 带 绿洲 
区 ,是 新 疆土 地 开发 活跃 区 域 。20 世纪 90 FARR 
以 来 随 着 国家 西部 大 开发 战略 的 实施 ,新 疆 也 进入 
到 新 一 轮 的 土地 开垦 与 生态 建设 时 期 ,耕地 面积 迅 
y sk 7 7? , 古 尔 班 通 古 特 沙 漠 南 缘 及 西部 均 有 土 
地 开明 ,沙漠 西南 部 地 区 的 开垦 最 为 明显 。 其 二 , 实 
地 调查 发 现 ,沙漠 边缘 地 带 也 普遍 实施 了 绿化 建设 ， 
重建 和 新 建 了 荒漠 林 和 防 沙 林 ,从 遥感 影像 上 也 可 
以 看 出 沙漠 南 缘 EVI 值 较 大 。 其 三 ,MOD13Q1 产 
品 的 空间 分 辩 率 为 250 m ,边缘 地 区 包含 绿地 的 混 
合 因 素 也 对 EVI 值 有 一 定 影响 。 

古 尔 班 通 古 特 沙漠 是 我 国 独 具 景 观 特色 的 固定 
半 固 定 沙漠 ,生物 .生态 系统 丰富 而 多 样 。 近 年 来 ， 
各 地 林业 .环保 等 部 门 不 断 加 强 对 其 管 护 ,许多 地 方 
建 起 了 围栏 并 设置 了 管 护 站 点 。 除 了 进行 一 些 必 要 
的 科学 研究 外 ,已 对 整个 沙漠 实施 了 全 面 的 保护 措 
施 ,这 一 系列 的 封 , 管 育 相 结合 的 综合 治理 措施 都 
促进 了 整个 沙漠 地 区 的 植被 保育 。 

古 尔 班 通 古 特 沙漠 EVI 的 年 内 变化 曲线 呈现 
先 上 升 后 下 降 的 趋势 。 从 3 月 开始 ,冬季 的 积 雪 开 
台 融 化 , 雪 水 消融 下 渗 , 土壤 得 到 了 水 分 的 补给 ,为 
植物 种 子 的 萌发 尤其 是 短命 植物 的 生长 发 育 创造 了 
条 件 ,4 月 初 短命 植物 开始 生长 ,并 在 接 下 来 的 2 个 
月 时 间 里 ,完成 整个 营养 周期 ,6 月 大 部 分 短命 植物 
干枯 死亡 ,此 时 , 梭 梭 ` 沙 拐 更 . 葛 类 等 植物 的 种 子 由 
于 春 冬 雨水 的 滋润 也 相继 萌发 生长 ,到 7 月 达到 一 
年 中 植物 生长 最 旺盛 的 时 间 ,随后 整个 沙漠 的 植被 
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5 结论 


(1) 2000—2018 年 , 古 尔 班 通 古 特 沙漠 整体 及 
其 不 同 沙漠 类 型 EVI 均 呈 显著 增加 的 趋势 ,这 表明 
近 19 a 上 古 尔 班 通 古 特 沙漠 的 植被 状况 向 好 的 趋势 
发 展 , 整 个 古 尔 班 通 古 特 沙漠 EVI 的 增长 速率 为 
0.016 0 - (103) (P<0.01)。 这 与 该 区 域 降水 、 
温度 等 气候 变化 以 及 人 工 生 态 建 设 和 加 强生 态 环 境 
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管 护 等 有 关 。 

(2) 生长 季 内 古 尔 班 通 古 特 沙漠 及 其 不 同 沙 漠 
类 型 的 EVI 变化 趋势 大 致 相同 ,EVI 从 3 月 开始 逐 
渐 增 大 ,并 在 6 一 7 月 达到 一 年 中 的 最 高 值 。 

(3) 古 尔 班 通 古 特 沙 漠 大 部 分 区 域 (96. 46% ) 
的 EVI S45 AU S, DURS 3.5496 KIRAY EVI 呈 减 
小 趋势 ,其 中 有 4396 的 区 域 显 著 增 大 (P <0.1)， 
0. 0696 的 区 域 显 著 减 小 (P<0.1),EVI 增 大 趋势 明 
显 的 区 域 主要 集中 在 古 尔 班 通 古 特 沙 席 的 南 缘 , 增 
大 趋势 最 为 明显 的 区 域 在 沙漠 西 部 边缘 农业 开 藤 区 
域 ,EVI 减 小 的 区 域 主 要 分 布 在 沙漠 北部 。 

(4) 古 尔 班 通 古 特 沙 漠 EVI 波动 较 大 的 区 域 
集中 在 沙漠 西 缘 与 南 缘 , 其 中 波动 最 大 的 区 域 在 沙 
漠 西 部 边缘 靠近 克拉 玛 依 的 开荒 区 域 , 沙漠 腹地 
EVI 波动 程度 较 小 。 

(5) EVI 指数 能 够 较 好 地 反映 以 固定 和 半 固 定 
为 主 的 沙漠 区 的 植被 覆盖 变化 ,其 反映 的 植被 状况 对 
于 沙漠 地 貌 类 型 的 空间 划分 具有 重要 的 参考 意义 。 
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Spatiotemporal Change of EVI in the Gurbantunggut Desert from 2000 to 2018 
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Abstract: Based on the MODIS EVI data,the spatiotemporal change of EVI in the Gurbantunggut Desert was 
studied using the one-dimensional linear regression analysis , empirical orthogonal function ( EOF) and variable co- 
efficient. The results showed that the interannual variation of EVI in the Gurbantunggut Desert from 2000 to 2018 
was in a significant increase trend with an increase rate of 0.016 0 - (10a) "' (P «0. 01) ,and it on fixed dunes, 
semi-fixed dunes and mobile dunes was also in a significant increase trend. The variation trend of EVI in the whole 
desert and different types of the deserts was similar in growth season. From March, EVI increased gradually and 
reached the highest value during the period from June to July. The EVI was in an increase in most areas of the Gur- 
bantunggut Desert , especially in the southern and western marginal zones,the increase trend was low in the hinter- 
land of the desert , and the areas with a decrease of EVI were mainly distributed in the northern part of the desert. 
The areas with high EVI fluctuation were mainly distributed in the southern and western marginal zones of the Gur- 
bantunggut Desert ,the fluctuation of EVI was the most significant in the western marginal zone , and it was low in the 
hinterland. EVI can used to reflect the change of vegetation coverage in the desert with fixed and semi-fixed dunes, 
and these results have important reference significance for spatial division of regional desert types. 
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